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We factorize the infrared divergences of the non-forward BFKL Green function for a general 
f-channel projection of the color quantum numbers and study the properties of the infrared finite 
remainder in the case of color octet exchange. The octet Green function is compared with the singlet 
case for different values of the momentum transfer. The octet Green function plays an important 
role in the calculation of the finite remainder of the two-loop six-point MHV planar amplitude in 
M = ^ SYM as it was demonstrated by Bartels, Lipatov and one of us in 1 1 , 2|. A comparison with 
regularizations preserving conformal invariance at large momentum transfer is shown. 

1 Introduction 

The Balitsky-Fadin-Kuraev-Lipatov (BFKL) formalism |I3}I3 has been used to investigate the prop- 
erties of scattering amplitudes in Quantum Chromodynamics (QCD) and M = A supersymmetric 
Yang-Mills theory (MSYM) in certain kinematic regions (multi-Regge (MRK) and quasi-multi- 
Regge kinematics) where logarithms of the center-of-mass energy are enhanced. In the leading and 
next-to-leading logarithmic EH approximations the elastic 2^2 and inelastic 2 2 + n am- 
plitudes have an iterative structure given in terms of reggeized gluon propagators and "squared" 
effective emission vertices. This iterative structure can be written for a general color representation 
in terms of a gluon Green function obeying a linear Bethe-Salpeter equation. The so-called "hard" or 
BFKL pomeron corresponds to the projection in the <-channel of a singlet in color space. This is the 
most widely used representation since it is related to physical observables such as hadron structure 
functions at small values of Bjorken x in deep inelastic scattering or inclusive dijet production with 
a significant rapidity separation at the Large Hadron Collider ll8] 4T6l . 

The BFKL resummation program has also been very useful from a more formal point of view. It 
was in a generalized leading logarithmic approximation, iterating the BFKL kernel in the s-channel, 
where the Bartels-Kwiecinski-Praszalowicz (BKP) equation was proposed |17 18 1 and found to 
have a hidden integrability |fT9l42T1l . being equivalent to a periodic spin chain of a XXX Heisenberg 
ferromagnet ll22l[24l . This was the first example of the existence of integrable systems in QCD. A 
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similar integrable spin chain, an open one this time, has recently been found by Lipatov in kinemat- 
ical regions of n-point maximally helicity violating (MHV) and planar {Nc oo) amplitudes in 
MSYM where Mandelstam cut contributions are maximally enhanced 1125 11261 . Corrections to the 
Bern-Dixon-Smirnov (BDS) iterative ansatz lIZTl for this class of amplitudes were found in MRK 
by Bartels, Lipatov and one of the authors in the six-point amplitude at two loops in [T,'2l. These 
corrections have been understood as part of the finite remainder [28-ri30J to the amplitude which 
corresponds to the anomalous contribution of a conformal Ward identity. A very compact formula 
for the amplitude can be found in terms of Goncharov polylogarithms 131113211 which is in agree- 
ment with the corresponding Regge limits 1, 3 3,-/3 5 J (also at three loop level ll36l ). Very recently, 
Romagnoni and one of us have found a novel relation between the BFKL equation in the forward 
limit and the sl(2) invariant XXX spin -1/2 chain |3T|. 

The octet gluon Green function is a fundamental ingredient of the finite remainder of scattering 
amplitudes with arbitrary number of external legs and internal loops. We investigate its all-orders 
structure in detail in the present work. In Section 2 we present an iterative representation for the 
Green function with a general color projection in the t-channel written in rapidity and transverse 
momentum space. We show how its infrared divergencies can be factorized in a simple form, leaving 
an iterative finite remainder that we investigate in detail in Section 3. Here we use advanced Monte 
Carlo techniques [38 J to study the exclusive information present in this Green function in both the 
singlet and octet projections. There appears faster convergence in the octet solution with a similar 
behavior in both forward and non-forward solutions. The collinear regions, related to the calculation 
of anomalous dimensions, are also investigated. In Section 4 we discuss the factorization of the 
infrared divergencies used by Bartels, Lipatov and one of us in Ref. [!l , 2] where they are extracted 
introducing a momentum transfer dependent gluon Regge trajectory. A similar regularization is used 
by Fadin and Lipatov in Ref. 1391 , which is also discussed here. We show the relation among these 
representations and discuss the change in the role played by the n = and n = 1 conformal spins. 
In Section 5 we present our Conclusions and scope for future work. 



2 The Green function for a general colour projection 

The non-forward BFKL equation for a general projection of the colour quantum numbers in the 
t-channel can be found, at leading and next-to-leading order, in Ref. II391I40II . In the present work 
we are interested in comparing the gluon Green function in the singlet case with that of the octet 
projection. The only difference between both solutions is in the "real emission" part of the kernel, 
which in the color octet case carries an extra factor of 1/2 with respect to the singlet case. In the 
octet case this spoils the complete cancellation of infrared divergencies present in the singlet, or 
pomeron, projection. 

To show that the extra infrared divergencies that appear in the non-singlet representations can be 
written as a simple overall factor in the gluon Green function, we regularize half of the divergencies 
in the gluon Regge trajectory using dimensional D — i — 2 e regularization, while the remaining 
ones are treated using a mass parameter A, which is also used to regularize the phase space integral 
of the "real emission" sector. The dependence on A will cancel out while the dependence on e will 
remain in the factorized term. We use the notation q'j = — q, where q is the momentum transfer 
and all two-dimensional vectors are represented in bold. 

In this way, the non-forward BFKL equation for a kernel in a general color group representation 
TZ in the t-channel can be written as 
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as 



A2 



|a^(qi,q2;q) = (5(2) (q^ - q^) 
qi2(qi+k)2-q2k2 



(q'l + k)2q2 



t;,^(qi +k,q2;q). (1) 



Note that compared to |40| all we have done is to add and subtract two gluon trajectories with a c-jz 
factor in front of them. The solution of this equation is independent of A in the A — limit, /i is 
the renormalization scale in the MS scheme. As we pointed out before, in the singlet representation 
ci = 1 and all infrared divergencies cancel out. In the octet representations cg^ — cs, =1/2 and the 
divergent behaviour of the gluon trajectory times {c-ji — 1) remains. Other representations, which 
we will not investigate here, are cio = Cjg = and C27 = — CAr^>3 = — l/(4A^c)- 
We solve Eq. ([T]) first introducing the following notation: 



{ 



UJ — UJ 



Hqiiq)}^;^ (qi;q2;q) = '5(^)(qi-q2) 

+CTZ J ^^(k2-A')e(qi,k; q) (qi + k, qa; q) . 
This can now be iterated and, performing the Mellin transform, 

/dcu 
^e"^"^^,^ (qi,q2;q) , 

we finally obtain 

J-(qi,q2;q;Y) = cxp {c^(^^^)(qi; q)Y} 1^(2) (q^ _ q^) 

00 n „ ,21 / "^"1 \ / ^ 

+ En^^ / ^e(k.^-A^)dqi + Ek^k.;qM''Mqi + Ek^-q2 

n=li=l J ''^ I \ 1 = 1 / \ 1=1 / 



(2) 



(3) 



dyi exp • 



UJ 



qi + X!'''!q 



qi +X!^''q 



1=1 



1=1 



(4) 



where the rapidity yo = Y ^ log s is given by the logarithm of the center of mass energy -/i. 

The A and 1/e dependence in the Reggeon propagators (i.e., all the exponentials in Eq. (|4]i) 
appear only in the initial multiplicative factor since 



qi 



1=1 



as 
-ylog 



(qi + EUk^)'(qi + EUki 
(qi + E;iik,)'(q'i + E:iik; 



(5) 



for any t-channel colour representation. 

This implies that the gluon Green function can then be written in the form 



-^(qi,q2;q; Y) 



Vq?qf 



-a,Y 



^2(1-Ctc);^2ctc 



,5(2) (qi - q2) 
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+En 



)'(qi + E;;!k()^ 



rVi-i 
X / dyi 
Jo 



(qi + Eii!k.)'(qi + Eii!ki 



2\ -tVi 



(qi + ELik;) (q'l + EUk^)^ 
The forward hmit, q = 0, is particularly simple: 



qi 



q2 



(6) 



V 



•^(qi,q2;Y) = e 



- p(l-C7i)^Y 



q? 



^2(1-ctc);^2ctc 



-a,Y 



<5(^^(qi-q2) + En^«-^/^^(k?- 



fVi-l 

Jo 



qi + ELl ki 



2\ O'sVi 



V (qi + El=i ki) 
In the case of colour singlet ci = 1 the solution is well-known: 

dui f d'-f f qf 



(7) 



1=1 



27ri Vqi/ w-asXn(7)' 



where 



X„(7) = 2*(1)-* 7 + ^ -^1-7+^ 



(8) 



(9) 



is the eigenvalue of the BFKL kernel and ^' is the logarithmic derivative of Euler's Gamma function. 
6 is the azimuthal angle between the two-dimensional vectors qi and q2. 

From now on we wiU focus on the description of the infrared finite remainder in the form 



^(qi,q2;q;Y) = J" (q^, q^; q; Y) |^-^=^ j 



aa(cTC-i)Y 



= (- 



<5('' (qi - q2) 



J n (qi+ELi 



(qi+ELikO^-q^k2^ 



n=l i=l 



rvi-i 
X / dyi 



^i? (qi+Eiliki)^ 



(qi + Eii!ki)'(q', + E;i!k,) 

(qi + Elik;)'(q'i+Elik, 



2\ -^Vi 

^^^^ (qi+f^k,-q2 ) I, (10) 



1=1 



where we have removed the factor with the e and dependence. We have checked that this func- 
tion "H is A independent for small A. In the following section we will compare this infrared finite 
remainder H for ctj = 1 and 1/2. 
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3 Comparing the octet with the singlet gluon Green function 

An interesting question is to study the convergence of the sum defining the function T-l in Eq. ( fTOb . 
For a fixed value of Y and the coupling as a finite number of terms in the sum is needed to reach a 
good accuracy for the gluon Green function. As the value of the effective parameter OgY gets larger 
the Green function is more sensitive to high multiplicity terms, following a Poissonian distribution as 
it can be seen in Fig.II] In both plots at the top of the figure, which correspond to the singlet forward 
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Figure 1 : Distribution in the contributions to the BFKL gluon Green function with a fixed number 
of iterations of the kernel, plotted for different values of the center-of-mass energy, and a fixed 
as = 0.2. 

q = (left) and non-forward, with q = 5 GeV (right), cases, we can see how the distribution in the 
number of iterations of the kernel gets broader for larger center-of-mass energies. Despite this there 
is good convergence, e.g., for Y = 4, gi = 15 GeV and q2 = 30 GeV, it is enough to keep up to 16 
terms in the sum to get a very accurate result. It is remarkable that the convergence is much better in 
the octet case where it is possible to get the Green function with a small number of terms, see Fig.[T] 
(bottom). This can be qualitatively understood if we think of the "real emission" terms as pushing 
the Green function to grow with Y. In the case of the singlet they carry the same coefficient as the 
contributions from the gluon Regge trajectories, which make the Green function to decrease with Y. 
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For the octet this balance is broken in favor of the latter making the Green function to grow much 
slower with an effective reduction in the number of emissions. 




Figure 2: Projection of the gluon Green function on different Fourier components in the azimuthal 
angle between the qi and q2 transverse momenta. 

It is also instructive to compare both solutions in terms of the different Fourier components 
in the azimuthal angle between the two momenta qi and q2. In the particular case of the singlet 
and forward case this expansion is given in Eq. A complete analysis for the singlet and octet 
projections is shown in Fig. |2] There are no surprises since both solutions behave quite similarly, 
with the only rising component being the n = one. It is possible to appreciate the similarity 
between both color projections in the full range of azimuthal angles in Fig. |3] 

The most striking difference between the singlet and octet cases appears in the collinear limit. 
This is investigated in Fig. |4] In the region with low transverse momenta the octet Green function 
gets a big decrease when compared to the singlet one. This indicates that the corresponding anoma- 
lous dimensions governing this limit are quite different in both cases. The collinear behaviour for 
the octet representation at large momentum transverse is shown in Fig.|5] 
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Figure 3: Gluon Green function dependence for the full range in azimuthal angles. 
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Figure 4: ColUnear behaviour of the gluon Green function. 
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Figure 5: Collinear behaviour of the octet gluon Green function at large momentum transfer 
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4 Conformally invariant representations 

As it was shown in Li, 2J, in order to have an expHcitly SL{2, C) invariant representation of the 
gluon Green function at large momentum transfer it is possible to subtract in Eq. ([T]) the gluon 
Regge trajectory depending on the total momentum transfer q^. To translate that representation to 
our present notation we need to work with the following modified Green function J: 



J'(qi,q2;q;Y) 



2 \ 2cTi 



? /2 

qiQi 



<^('' (qi-q2) 



A^TT^_ / ^ ^*fln,2 \2^ 



n— 1 i—1 



qi+E/=ik/ (qi+ELikO'-q'kf 



(q'l + EUkOMqi + E/iik, 



Vi-^- I (qi+E/=!k,) (qi+E/=!k/ 

^qi + ELik,)'(qi + ELik; 



S(2) 



qi 



1=1 



I - q2 



(11) 



The treatment of the infrared divergencies is equivalent in Eq. ( fTOl i and Eq. ( fTTT i since they carry 

the same multiplicative factor A^'^^"^^. They are different in a relative normalization of the form 

N (l-c7j.)asY 

/q^q^/ q^ J which is not relevant in the singlet case but diverges when q^ — for the 



octet representation. In the "real emission" part the modification with respect to Eq. ([T]i is 



1 + 



q'i2(qi+k)2-q2k2 



1 



qf(qi+k)2-q2k2 



q?qf 



(12) 



(q'i+k)2q2 ^ (q'i+k)2q2 ^ (q^ + k)2(qi + k)^ ' 

In the numerical analysis of the solution to this representation of the BFKL equation we find agree- 
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0.01 " 




Figure 6: Left: Projection of the function J' on different conformal spins. Right: Collinear be- 
haviour of the function J' for two values of the rapidity Y. 

ment with the work in Uill where it was shown that at large q^ the dominant conformal spin is 
n = 1, since it is the only one having a positive asymptotic intercept. The intercepts for all the other 
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conformal spins, including n = 0, are negative. We show this behaviour in Fig. |6] (left) where we 
find that, in the limit ~ q^, the n — 1 projection in the azimuthal angle is the only one 

growing with energy. 

The collinear structure of the solution in the form of the function J' is illustrated in Fig.|6](right) 
for two different values of Y. We find agreement with the one obtained in the previous section. It is 
also noteworthy the overall decrease of the Green function as Y grows which occurs because in this 
plot we show the azimuthal angle integrated solution for which only the n — component survives. 

In a recent work Fadin and Lipatov [39] have presented the NLO non-forward BFKL equation in 
the adjoint representation for the JV = 4 supersymmetric theory. Their normalization for the Green 
function and treatment of the infrared divergencies, showing for simplicity only the LO pieces, is 
the following: 



a 



a, f d^k 

— ^ / ( 

4 / 7rk2 



(k2 - 



qf(qi+k)2-q2k2 
(q'l + k)2q2 



(qi + k) 



r^?^(qi+k,q2;q). (13) 



The corresponding solution in rapidity space reads: 



'B(qi,q2;q; Y) 



q2A2^ 

q?qf . 



#Y 



5^'^ (qi - q2) 



n=l i—1 



qi + ElJi k; 



qi + ELi k) 



1+^ 



q'l+EUk, (qi+EUkO^-q^k? 



(q'i + ELikOMqi + E;;ik, 



' (q; + : 

(q; + EU: 



dyi 




The difference with the J' function is the normalization of the Green function, while the treatment 
of the infrared divergences is the same. It is then natural to obtain a very similar behaviour for B as 
it is shown in Fig.|2l 



5 Conclusions and scope 

In this short letter, a detailed study of the solution to the octet BFKL equation has been carried 
out. So far, this solution has only been discussed in ||Tl|2l[39l where it naturally appears in the 
context of maximally helicity violating planar amplitudes in Af — A super Yang-Mills theory. In 
certain kinematical regions, the gluon Green function in the octet sector plays a fundamental role 
to obtain the so-called "finite remainder function". It is therefore very important to understand its 
structure in depth, not only to the first orders in the 't Hooft coupling, but also its all-orders features. 
This has been the target of the work presented here. We have discussed several representations 
directly written in the transverse momentum space which allow for the factorization of its infrared 
divergencies. The dominant Fourier component in the azimuthal angle is different in these 
representations, also showing a different collinear limit when compared to the singlet case. The 
integration of these representations with the corresponding impact factors to investigate complete 
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Figure?: Left: Projection ofthe function Son different conformal spins. Right: Collinear behaviour 
of the function B for two values of the rapidity Y. 

scattering amplitudes to all orders in the coupling will be the subject of our future investigations. 
The study of the next-to-leading order color octet kernel is also underway. 
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